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Abstract 



We examine a nonlocal interaction that results from expressing the QCD 
Hamiltonian entirely in terms of gauge-invariant quark and gluon fields. The 
interaction couples one quark color-charge density to another, much as elec- 
tric charge densities are coupled to each other by the Coulomb interaction in 
QED. In QCD, this nonlocal interaction also couples quark color-charge densi- 
ties to gluonic color. We show how the leading part of the interaction between 
quark color-charge densities vanishes when the participating quarks are in a 
color singlet configuration, and that, for singlet configurations, the residual 
interaction weakens as the size of a packet of quarks shrinks. Because of this 
effect, color-singlet packets of quarks should experience final state interactions 
that increase in strength as these packets expand in size. For the case of an 
SU(2) model of QCD based on the ansatz that the gauge-invariant gauge field 
is a hedgehog configuration, we show how the infinite series that represents 
the nonlocal interaction between quark color-charge densities can be evalu- 
ated nonperturbatively, without expanding it term-by-term. We discuss the 
implications of this model for QCD with SU(3) color and a gauge-invariant 
gauge field determined by QCD dynamics. 
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I. INTRODUCTION 



In previous work, [I|-|3| we developed a gauge-invariant formulation of QCD in which the 
QCD Hamiltonian is described as a functional of gauge-invariant operator-valued quantities 
- spinor (quark) fields, gauge fields and their canonical momenta. One outgrowth of this 
work was the demonstration of an important similarity between QED and QCD. In a for- 
mulation of QED in which the charged spinor field represents a gauge-invariant quantity, 
a nonlocal interaction between charge densities — the Coulomb interaction — appears as 
part of the Hamiltonian. On the other hand, as is well-known, there is no explicit nonlocal 
interaction in standard versions of QED in covariant gauges in which the charged spinor, 
ip, represents a gauge-dependent field. However, when covariant-gauge QED is transformed 
to a representation in which Gauss' law is implemented and ip becomes a gauge-invariant 
charged field, the same nonlocal Coulomb interaction that is seen in the Coulomb gauge ap- 
pears explicitly in that case as well, even though the covariant-gauge condition continues to 
apply to it. |4[] These observations illustrate that, when only gauge- invariant fields are used 
in constructing the Hamiltonian — in a variety of gauges, including but not limited to the 
Coulomb gauge — the interactions between charged fields and pure-gauge parts of gauge 
fields vanish, and a nonlocal interaction between charge densities appears in their stead. 
The fact that the Coulomb interaction is, by far, the most important electrodynamic force 
in the low-energy regime, provides strong incentive for formulating QCD in terms of gauge- 
invariant fields, to explore the implications of these observations for our understanding of 
QCD. 

The organization of this paper is as follows: In the section following this Introduction, we 
review and expand on relevant material from our previous work, which provides a foundation 
for describing QCD in terms of gauge invariant fields. This part of our discussion is exact, 
and requires no approximations. In a later section, we will use these results to suggest 
how quark confinement and color transparency could be understood as a consequence of the 
gauge-invariant formulation of the QCD Hamiltonian. And finally, we will invoke an ansatz 
about the form of the gauge-invariant gauge field to illustrate a technique for evaluating the 
nonlocal interaction nonperturbatively. 



II. THE QCD HAMILTONIAN AS A FUNCTIONAL OF GAUGE-INVARIANT 

FIELDS 

We will review here some technical developments that we discussed previously, and 
that are essential for the investigation reported in this work. One of these developments is 
the construction of a set of gauge-invariant quark and gluon operator- valued fields. Another 
is the transformation of the QCD Hamiltonian to a representation in which it is expressed 
in terms of these gauge-invariant fields. In our work, the QCD Hamiltonian, H, is expressed 
in terms of gauge-invariant operator-valued fields in a representation in which ijj designates 
the gauge-invariant quark field. It is given by 
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H = dr 



±n?(r)II«(r) + i^.(r)^.(r) + ^(r) (p m _ l(Xidi) ^( r ) + . ( 2 .1) 



if' describes interactions involving the gauge-invariant quark field, and can be expressed as 

H' = H J . A + H g + H LR . (2.2) 

As shown in Ref. [|J, Hg vanishes in the representation (the so-called M representation) in 
which this Hamiltonian is described, and will not be given any further consideration in this 
work. Hj-A describes the interaction of the gauge-invariant gauge field with the transverse 
gauge-invariant quark color-current density, and is given by 



Hj-A = ~9 J dv^(r)a4^(r)A h GH (r) 



(2.3) 



Hlr is the nonlocal interaction 



Hlr = Ha-Q + Hr 



(2.4) 



withfl 



Hg-C 



and 



drl +Tr 



J29 r+1 (-lYflff^Ut(r) f Al H {v) VZ\t)%V c {t)& (^(r)^(r)) 



r=0 



+Tr 



r=0 



(2.5) 



%Q- IEE 9 r+r '(-ir r 'f^ff^'§ (^totfto) § (s&totf'to) ; (2.6) 



r=0r'=0 



here Ag U (r) is the transverse, gauge- invariant gauge field constructed in Ref. || and j£ is the 
gauge-invariant quark color-charge density j'q = gip Ji (\ h /2)tp ; /ff 7 is the chain of structure 
constants 



_ j a [l]/36[l] j6[l]a[2]b[2] jb[2]a[3]6[3] . . . j6[r?-2]a[r?-l]6[r;-l] yftfa-ljafafr 



(2.7) 



summed over repeated Lie group indices; V c (r) -yVc(r) is a quantity that transforms like an 

SU(N) vector (where iV = 3 for QCD and N = 2 for Yang-Mills theory); and (T^(r)jft(r] 
is given by 



^f)(r)jo(r) = 4« (1) (r) ^ (• • • (f Q (r))))) . (2.8) 



Eq. ( |2.6| ) can be understood as the non-Abelian analog of the Coulomb interaction in QED, 
and to have a structure similar to that of the Coulomb interaction, when the quantity 



x Eq. ( [2.5D corrects a typographical error in Eq. (23) in Ref. g. 
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^Q)( r ) = Sr=o ffr) h (~9) r (^(f)( r )io ( r )) IS substituted for the Abelian j (r) that represents the 
electric charge density in QED. The same /Q[(r) also participates in the nonlocal interaction 
described in Eq. (|2.5|) , where it couples to "glue" -color, given by 

i4& 4 (r)Il*(r). (2.9) 

Before proceeding with an application of our earlier work to a discussion of quark color 
confinement and color transparency, we need to clarify some questions about the formalism 
we have constructed. Our first remark addresses the fact that gauge-invariant fields are not 
generally unique, either in Abelian theories, such as QED, or in QCD. Since the Gauss' law 
operator is the generator of infinitesimal gauge transformations, any unitary operator that 
commutes with the Gauss' law operator can be applied to a gauge-invariant state, or can be 
used to unitarily transform gauge-invariant operators, without interfering with their gauge 
invariance. A number of such unitary operators that commute with the Gauss' law operator 
for QED were displayed in Ref. ||], and further operators, that transform gauge-invariant 
operators in other Abelian gauge theories without interfering with their gauge invariance, 
have also been constructed, p Constructing such operators for QCD would be more difficult, 
but there is little doubt that it would be possible to do so. We therefore need to address 
the non-uniqueness of the gauge-invariant operators used in our formulation. 

In previous work, HH, we have used the following criterion for accepting a gauge- 
invariant field as useful for a gauge-invariant formulation of a gauge theory; we propose 
to apply that same criterion in this case as well. We require that when the interactions 
between the pure gauge degrees of freedom and gauge-invariant matter fields have been 
eliminated in favor of nonlocal interactions involving gauge-invariant matter fields — Hq^q 
and H g _Q in the present case — the remaining interactions are restricted to interactions of 
the gauge-invariant gauge field with transverse current densities only. The unavailability of 
a longitudinal component of the gauge-invariant gauge field in Hj-A precludes the forma- 
tion, through the operation of current conservation, of a coupling between charge density 
and a longitudinal gauge field. No further interactions between matter-field charge densities 
can therefore be transmitted through virtual loops of gauge-invariant gauge field compo- 
nents, making the nonlocal interactions the dominant features in a description of low-energy 
dynamics. 

In most gauge-invariant formulations of QED, such as those in the Coulomb gauge, the 
Coulomb interaction, — |/c/rjo(r)V _2 jo( r ), and the interaction of the transverse (gauge- 
invariant) part of the gauge field with the current density, constitute the interaction Hamil- 
tonian. That formulation satisfies the criterion we have proposed. In axial (A 3 = 0) gauge 
QED, however, the nonlocal interaction is given by — ^Jdrj (r)d 3 2 jo(r) and is accompanied 
by further interactions of the gauge field with charge as well as current densities. Further 
unitary transformations are then necessary to eliminate contributions from these additional 
interactions, through virtual photon loops, to the formation of forces between static charges. 
Appropriately chosen unitary transformations not only eliminate these interactions among 
static charges through virtual loops, but also transform — |/ drj (r)d^ 2 jo(r) to the Coulomb 
interaction — |/drjo(r)V~ 2 jo( r ) and thereby restore rotational symmetry to the nonlocal 
interaction, and still maintain the axial gauge condition. |4[] It should also be noted that, 



Kj(r) 



gf dae Jr 



-V c (r 
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in QED, our proposed criterion requires that gauge- invariance be imposed on matter fields 
through unitary transformations that do not dress charged fields with transverse propagat- 
ing photons that have no role in implementing Gauss' law. Q A similar requirement, that 
charged particles not be dressed spuriously with transverse, propagating photons, was pro- 
posed by Haagensen and Johnson as necessary for avoiding "false confinement" in Wilson 
loop calculation for QED. || 

When applied to QCD, we note that the nonlocal interactions in H' are Hq_q and 
Hg_Q, and that Hj-A describes the interaction of the transverse, gauge-invariant gluon field 
^gi j( r ) w ith quark color transverse current densities only. Hg is a term in the transformed 
Hamiltonian in which Q a appears either on the extreme left or right. Q a represents the 
Gauss's law operator in the representation in which the transformed Hamiltonian H is 
expressed. Q a always annihilates either the "bra" or "ket" state vector of matrix elements 
taken between allowed states. And, since Q a commutes with H, state vectors that implement 
Gauss's law initially, will continue to do so as they time-evolve under the influence of H. 
The nonlocal Hq_q and H g _q therefore describe all the interactions between static quarks. 
No interactions between quark color charge densities can be transmitted through virtual 
loops generated by the gauge- invariant gluon field Aq H (y) without involving Hq_q or H 9 -q; 
and, therefore, Hj-A can be expected to make only relatively unimportant contributions 
to low-energy processes. This fact supports our choice of gauge-invariant quark and gluon 
fields as appropriate for formulating a gauge-invariant description of QCD dynamics. 

We will also make some clarifying remarks about the relation between two representations 
used in this and earlier work: the C and the M representations. As we pointed out in 
earlier work, 0]|| the Gauss's law operator Q a {r) = diH^(r) + gf abc A b (r)U^(v) + jo( r ) with 
j'o( r ) = 9 V^( r ) 4f ^(r) , and the "pure glue" Gauss's law operator Q a (r) = <9jIT/(r) + 
gf abc A b (r)U^(r) are unitarily equivalent, so that Q a may be taken to represent Q a (r) in 
a different, unitarily equivalent representation. We refer to the representation in which 
Q a (r) is the Gauss's law operator, the C representation; and the representation in which 
Q represents the entire Gauss's law operator, with the color-charge density Jq(t) included 
(though implicitly only) the M representation. The unitary equivalence is expressed as |J 

g a {v)=u c g a {v)u c l . (2.io) 

The Gauss's law operator is the generator of infinitesimal gauge transformations; and the 
criterion for the gauge invariance of an operator £ is that it commute with the Gauss's law 
operator. In this work, when we use the Gauss's law operator to determine gauge invariance, 
it is important to distinguish between the Gauss's law operators in the C and the M repre- 
sentations. An operator £ represents two different quantities in the two representations. The 
quantity represented by £ in the C representation has the form U^ 1 ^Uc in the J\f represen- 
tation. But £, when appearing in the M representation, refers to a different quantity, whose 
form in the C representation would be UcC^c ■ Since the quark field ip trivially commutes 
with G a (r), ip is manifestly gauge- invariant in the M representation. The unitary operator 
Uc transforms the quark field ip so that 

Vci(r) = Uci>{r)Uc X = V c (r)1>(r) (2.11) 
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is the form that the gauge-invariant gauge field takes in the more usual C representation. 
The unitary transformation given in Eq. (|2.11 ) has no effect on the gauge field. But the 



transformation it effects on the quark field is significant, because it allows us to use the quark 
field ip{r) to represent the gauge-invariant quark field in the TV representation. The unitary 
operator Vc(r) is a functional of gauge fields and of the Gell-Mann matrices X h ; its structure 
was discussed extensively in Ref. @. Since, in the Af representation, ip(r) is the expression 
for that same gauge-invariant field that is described by ipG\(r) in the C representation, Vc(r) 
is implicitly included in the quark field ip in the M representation. In the M representation, 
ip(r) therefore implicitly consists of glue as well as quark field components. When the 
quark field is gauge-transformed within the C representation, ip(r)— >exp[ig(X h /2)x h {r)]ijj(r), 
where x h { r ) is an arbitrary time- independent c- number field in the adjoint representation 
of SU(3); that same gauge transformation transforms Vc(r)— »Vc(r)exp[— ig(\ h /2)x h (r)}, so 
that ipG\( r ) is gauge-invariant. In the jV representation, the gauge invariance of ^(r) is a 
trivial consequence of the structure of the Gauss's law operator. 

These considerations have important consequences for the behavior of charge and current 
densities under gauge transformations. As a simple illustrative example, we consider the 
electron field operator ip in QED in the temporal gauge. Under a gauge transformation, 
A^^A^ — d^Xi an d the corresponding change in ip is ip— >exp(iex)ip ■ There are a number of 
ways of constructing gauge-invariant electron field operators in QED, || but one that is very 
useful for this discussion was provided by Dirac, || who defined the gauge-invariant electron 
field ipG\ = exp[— ie(l/d 2 )diAi } ip. Under a gauge transformation, compensating changes 
occur in ip and in (l/d 2 )diAi, so that ipQ\ remains gauge-invariant. The electric charge 
density operator can be represented either as eipQ\4>G\ or as eifjty. Since the longitudinal 
photons (they are zero-norm ghosts) used to dress the electron field to make it gauge- 
invariant are electrically neutral, they do not affect the electric charge density, and the two 
expressions, eipQ { ipQ\ and eip^ip, are identical. 

To discuss the non-Abelian case, we will use the SU(2) version of QCD, because the 
difference between the Abelian and the non-Abelian cases can be illustrated very graphically 
in that system. We therefore use the SU(2) version of Vc(r) in Eq. ( |2.11 ), which is given by 



Vc(r)su{2) = exp[— ig{f /2)-Z{y)}, where Z(r) = (1/d )diB i} and where Bi is a complicated 
functional of gauge fields, all of which commute with each other and with the spinor field 
ip. We can identify the gauge-invariant color charge density in the C representation of 
this SU(2) model of QCD as [j ]gi — S^aC^/^) V'gi- m this non-Abelian case, [j ]gi is no 
longer identical to j . When the substitution ipQ\ = exp[— ig(r/2)-Z(r)] ijj is made, the 
gauge-invariant [j ]gi can be expressed in terms of j = gip\f/2) ip by 

[Jo]gi = (i-Jo) Z + Zx Q xZ) cos(Z) + Q xZ) sin(Z) , (2.12) 



where Z is the unit vector Z(j) = Z{j)jZ(j). Eq. (|2.12|) shows that in this case, and in its 



SU(3) version, dressing the gauge-dependent quark field with the gluons required to make 
it gauge-invariant, does affect the color charge density, since the gluons themselves carry 
color. Under the infinitesimal gauge transformation 5 At = diSx + gAjXSx, the change in 
jo is exactly compensated by the change in Z, 0,|7j so that [j ]gi remains untransformed, 
and maintains its orientation in the — in this instance SU(2) — color space. The fact that a 
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quantity is gauge-invariant does not mean that it has no preferred orientation in color space; 
it does mean that, whatever orientation it has, will not be altered by a gauge transformation. 
The same remark applies to the gauge- invariant gauge field Aq H (t) as well. 

In view of the relation between [j'o]gi and jo, it is particularly important to realize that 
jo is the form that [j'o]gi takes in the Af representation. The color-charge density jg(r) 
as well as the color current density jf (r) therefore are gauge-invariant operators in the 
Af representation, and implicitly include glue as well as quark ingredients to constitute 
gauge-invariant color-charge and color-current densities. It is important to emphasize this 
point. The quantity that is represented by j'q (r) in the C representation is different from the 
quantity that is represented by j'o(r) i n the Af representation. jo( r )> when it appears in the 
C representation, consists of quark field components only, and transforms gauge- covariantly 
(not gauge- invariantly) under an infinitesimal gauge transformation, in the form 

5 bo(r)] c _ rep = 9f abe [j»] c _ r Jx e ■ (2.13) 

But Jo(r), when it appears in the Af representation, is gauge invariant! A gauge-invariant 
field tensor can also be defined as0 

*gi«(*) = dMr) - «M a Gli (r) - 9 f abc A Gil (r)A Gij (v) (2.14) 

and then F G ^.(r)f = i^(r)V c (r)£ V c ~ \r) and F£ Hj (r)Fg Hj (r) = ^(r)^.(r) . Similarly, 
Kg(r) and n"(r)II"(r) are gauge-invariant quantities, so that all the terms appearing in H 
are gauge- invariant. In Ref. ||, we verified that the expressions we obtained for the gauge- 
invariant fields — gauge as well as quark fields — by expanding Vq to arbitrary orders, agree 
with the expressions arrived at perturbatively by Lavelle and McMullan. P,|10|1 

The fact that, in the Af representation, j'q (r) and jf (r) represent gauge-invariant charge 
and current densities respectively, enables us to use them to represent physical observables. 
This makes it very convenient to formulate this work in the Af representation. It would, in 
principle, have been possible to carry out the investigation we are presenting here in either 
the C or the Af representation. But the use of the Af representation is a very powerful 
tool for expressing the Hamiltonian in this gauge-invariant formulation. We will therefore 
proceed with this discussion, by using the Af representation to examine how the nonlocal 
interactions Hq_q and H g _q could provide a mechanism for the confinement of quarks 
when these are not in color-singlet configurations. A number of authors have discussed the 
importance of implementing Gauss's law in the quantization of QCD and Yang-Mills theory. 



HH16I And other authors have also speculated that Gauss's law might be responsible for 
color confinement in QCD, and that, therefore, gauge-invariant degrees of freedom might be 
necessary to make color confinement manifest. H6|)P, [10|,|17 -|22]| A somewhat similar expression 



to our Eq. (|2.6|) — later expressed as Eq. (|3.2|) — has been given by T. D. Lee P3[ and 



2 In this non-relativistic notation, V{ refers to a contravariant vector and di refers to the covariant 
derivative. This convention is used extensively throughout this work. 
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by Christ and Lee, |24]] but not with the same gauge-invariant gauge field A s Gi between the 
inverse Laplacians that appears in our work. The gauge-invariant gauge field Aq, positioned 
between the inverse Laplacians, is instrumental in providing for gauge invariance of the 



interaction Hamiltonian. In Ref. |24|], unitary operators (u and K) transform the spinor 
(quark) fields so as to bring them into compliance with Gauss's law in a formulation in 
which Gauss's law is assumed to already be implemented in the gauge sector. Our work, 
reported in Ref. 0, is based on the initial implementation of Gauss's law for the gauge 
sector, so that the functionals we construct, that correspond to the angle variables in u 
and U in Ref. |]24|| , are explicitly given in terms of operator- valued gauge fields. These 
explicit expressions, which we obtained by implementing Gauss's law for the gauge sector, 
are required to establish gauge invariance of the quark and gluon fields and of the color 
charge and current densities. In spite of the important differences between the non-Abelian 



analogs to the Coulomb interaction that appear in Refs. [pffi , |24|| and the one given in our Eqs. 
(j2.6|), (|3~2|)and ( |3.3| ), the fact that nonlocal interactions of the same general structure appear 
in these two treatments, which are quite different in their objectives and in the technical 
procedures used, supports the idea that such nonlocal interactions can have an important 
role in QCD dynamics. 

We do not give a proof, in this work, that the long-range interactions Hq_q and H 9 ^q 
confine quarks. We do present arguments, however, that these nonlocal interactions, which 
arise naturally when the QCD Hamiltonian is represented in terms of gauge-invariant fields, 
are interesting candidates for describing low-energy QCD dynamics, and that these nonlocal 
forces might have an important role in the confinement of quarks and the non-confinement 
of color-singlet configuration of quarks. 



III. QUARK CONFINEMENT, COLOR TRANSPARENCY, AND THE 

NONLOCAL INTERACTION 

We will here examine to what extent the structure of JCq(t) can be responsible for the 
confinement of quarks, or of wave packets composed of quarks. To facilitate this discussion, 
we will display JCq(t) in the expanded form 

KS(r)= -j'S(r) + (-,) f'^AlUm I ^1 x| Jo°M + 



4vr|r - Y(i)| 



+ (-g) n f S W bS W. ■ ./ S («-2) 5 (n-l) S (n-l)/ S (n-l) 5 (n) a ^ ( | 1 4 ?( r )5 i j | 

A « l iy(n - 3))de J 4vr|y (n _3)-y ( n- 2 )| AGI ^ 

% f 47Hv dy{n ~l — |4ri(y ( n-i)) d k f d * jg(x) + • • • . (3.i) 

J 47T]y( n _ 2 ) — y(n-i) | J 47r|y (n _i) - x| 



Su 

Q\j '{y(n-2)) X 



s 



Eq. ( p.8|) can be substituted into Eq. Q2.6|), and a partial integration carried out, so that 
(l/2)/dr§|/Cg(r)§iX:£(r) is transformed into -(1/2) fdr/C£(r) (1 / d 2 ) }C b (r) . 



The first question to consider is how the nonlocal interaction described in Eqs. ( |2.6| ) and 
( j3.1|) would depend on the color of the participating quark configurations. We turn our 
attention to Eq. (|2.6|) and represent it as 



Hq-q = ~\f dvlC b (r) (1) K*{t) = i/rfrdx Jo 6 (r)^(r, x) Jo a (x) . (3.2) 

jF ba (r, x) is a Green function that is the non-Abelian analog of 1/ (47r|r — x|), but which 
differs from its QED analog in that it does not only refer to spatial points, but also depends 
on the gauge-invariant gauge field A 5 GU (r). jF 6a (r,x) can easily be read from Eq. ( |3.1|) , and 
shown to be 



29 fn* /rrn^(y)^ 

J 47T r — V 47T V - 



47r|r — x| "' J 47r|r — y| 47r|y — x| 

J 4vr|r-y 1 | J 4vr|y 1 -y 2 | J 47r|y 2 - x| 



(-l) {n - 1) (n + l)g n f s m bs mf s m 5 w s w...f s («-v 5 M a / —^1 — -AXUy 1 )d i / , y2 r x 

J 4vr|r - yi| J 4vr|yi - y 2 | 



4,1 (y 2 ) dj f . | dys | ■ ■ • / --^ r 4u J (y») dirr 1 — r + • ■ ■ • (3-3) 

J J 47r|y 2 -y 3 | J 47r|y (n _ 1) - y n | 47r|y n - x| 



jF ba (r, x) can be seen to be symmetric under the combined interchange r x and a^b, 
when we assume that A s cu (y)—>0 when |y|— >oo. Eqs. (ft.2|) and ( |3.3| ) define a nonlocal 
interaction between gauge-invariant quarks or packets of quarks that is an analog of the 
Coulomb interaction in QED. J rba {r, x) is determined not only by the inverse Laplacians, 
but also by the matrix elements of the gauge- invariant gauge fields A 5 Q H (y) in the gluonic 
medium in which the quarks are immersed. Once the spatial dependence of jF 6a (r, x) in 
a particular state of quark-gluon matter is fixed, it describes a force that acts on quark 
color charges at points x and y respectively. Particularly in cases in which quark- ant iquark 
creation is supressed and in which the quarks are at rest, or nearly at rest — as in the case 
of heavy, static quarks — jF 6a (r, x) can play a very similar role in QCD to the Coulomb 
interaction in QED. 

We will apply Eq. (|3.2|) to a state consisting of two packets of quarks that are immersed 
in glue, but are well-separated from each other in the sense that the quarks in each of the two 
packets can be represented using a complete set of orbitals, and that the two sets of orbitals 
occupied by the quarks in the two packets have a negligible overlap. We will not attempt to 
make any quantitative models for how such separated packets arise, but will only explore the 
interactions of such separated packets once they have arisen. We will furthermore assume 
in this discussion that quark-antiquark pair creation and annihilation can be neglected or 
"quenched" . We consider the expectation value 

(qQ | Hq-q I Qq) = ( q k ■ ■ ■ qt Q n ■ ■ ■ Qt | H Q _ Q \q 1 ■■■ q k Qi ■ ■ ■ Q n ) , (3.4) 
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where q± • • • q k represents a set of quarks in orbitals Uj(r) and Q\ • • • Q n a set of quarks 
in orbitals £/,(r), so that both sets of orbitals are localized, and the overlap of Uj(r) 
and Ui(r) is negligible. In | q% ■ • ■ q k Q x • ■ • Q n ), a set of creation operators for quarks 
lii " " " Ik , Qi, ■ ■ ■ Qn , are applied to a state \g), which represents a gluonic medium that 
obeys the Gauss's law (in the Af representation) 

G a \g) = o. (3.5) 

In Ref. [Q, a procedure was developed for constructing states such as \g). Since, in the 
M representation, quark creation and annihilation operators commute with the Gauss's law 
operator Q a , multiquark Fock states can be constructed by applying quark creation operators 
to \g) without invalidating Gauss's law applied to the resulting multiquark Fock state. 

Because quark and glue field operators commute, and because of the negligible over- 
lap of of the quark orbitals Uj(r) and Ui(r), (qQ \ Hg_q\Qq) simplifies, through cluster 
decomposition, to 

(qQ | Hq-q \Qq)= J dvdxf ba (r,x)(q k - • - qi \fi(r)\q v • •?*><£„• • -Qx\j^)\Qv • -Qn) , (3.6) 

where jF 6a (r, x) represents the expectation value jF ba (r, x) = (g\ jF ba (r, x) \g). 

In evaluating J<ix jF ba (r, x)j'q(x), — one of the integrals required for Eq. ( |3.2| ) — we 
can expand about a point xo located where the orbitals Ui(r) are large in magnitude. The 
integral then can be expressed as the Taylor's series, in a kind of "color multipole" expansion, 

J rfx |^ a (r, x ) + XA^ir, x ) + ~ X, t X 3 did^r, x ) + ■ ■ •} j a (x) (3.7) 

where Xi = (x — xo)i and di = d/dxi . When we perform the integration in Eq. ( |3.7| ), the first 
term contributes jF ba (r, Xo) Q a , where Q a = J c/xj'q(x) — the integrated "color charge". 
Since the color charge is the generator of infinitesimal rotations in SU(3) space, it will 
annihilate any multiquark state vector in a singlet color configuration. Multiquark packets 
in a singlet color configuration therefore are immune to the leading term of the nonlocal 
Hq_q. Color-singlet configurations of quarks are only subject to the color multipole terms, 
which act as color analogs to the Van der Waals interaction. 

The scenario that this model suggests is that the leading term in Hq^q, namely 
<2 6 jF ba (r , x )Q a , for a quark color charge Q a at r and another quark color charge Q b 
at x , as well as H g _q, which describes the coupling between quark color charge and color- 
bearing gluonic matter, are responsible for the confinement of quarks and packets of quarks 
that are not in color-singlet configurations. In this scenario, the multipole terms, which 
are the only parts of Hq_q that affect the dynamics of color-singlet quark configurations, 
do not confine packets, but result in final-state interactions that act on color- singlets as 
these move through, or emerge from a gluonic medium. To confirm this scenario, it would 
be necessary to evaluate jF ba (r,x) and its spatial derivatives. This requires knowledge of 
the spatial dependence of Aq H (t). In Refs. [@,[§, relationships of A 5 GH (r) to other gauge- 
invariant (and gauge- covariant) quantities were established that would, in principle, allow 
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specification of A s Gii (r) (or, rather, the specification of the set of possible values of A GH (r), 
since the nonlinearity of the equations that determine those relationships may well signify 
that A GH (r) does not have a single unique value). But the explicit construction of A GH (r) 
is not within the scope of this current work. We will, rather, explore some features of the 
behavior of color-singlet quark packets as these move through a gluonic medium, assuming 
that jF fea (r, x) is a reasonably well-behaved function of r and x. 

In order to examine the dynamics of color-singlet quark configurations as these move 
through a medium consisting of gauge-invariant glue and quark matter, we will explore 
the effect of the multipole terms on "small" color-singlet packets. Small packets occupy a 
limited region of space in the interior of gluonic matter, so that jF ba (r, x) does not vary 
significantly over the spatial domain in which the packet functions Uj(r) and £/j(x) make 
sizable contributions. Given the assumption that jF fea (r,xo) varies only gradually within a 
volume occupied by quark packets, the effect of these color multipole forces on a packet 
of quarks in a color-singlet configuration becomes more significant as the packet increases 
in size. As small quark packets move through gluonic matter, they will experience only 
insignificant effects from the multipole contributions to Hq_q, since, as can be seen from 
Eq. (p.7|) , the factors X{ , XjXj, ■ ■ -, X^iy ■ X^ n ), keep the higher order multipole terms 
from making significant contributions to fdx.F ba (r, x)j'q (x) when they are integrated over 
small packets of quarks. As the size of the quark packets increases, the regions over which 
the multipoles are integrated also increases, and the effect of the multipole interactions on 
the color-singlet packets can become larger. This dependence on packet size of the final-state 
interactions experienced by color-singlet states — i.e. the increasing importance of final- 
state interactions as color-singlet packets grow in size — is in qualitative agreement with 
the characterizations of color transparency and color coherence given by Miller and by Jain, 
Pire and Ralston. p5| , pp| Eqs.Q3.2p, ( |3.3j ) and ( ^.7[ ) generalize the multipole expansion of a 



time-independent nonlocal interaction from the electromagnetic case, in which the leading 
term vanishes for a neutral object and leaves a residue of higher order multipoles, to QCD, 
in which color neutrality — an attribute of color-singlet states — corresponds to electrical 
neutrality in electrodynamics. In the language of Jain, Pire and Ralston, [^] our analysis 
provides a model for how color neutrality protects a quark packet from the color-monopole 
force, with the result that such color-singlet packets can survive — when the region over 
which the multipole interaction is integrated is small enough and the effect of jF ba (r, x) on it 
coherent enough — to become asymptotic states observed in studies of exclusive processes. 

As previously mentioned, quark and gluon color is also coupled directly, in the form 

H 9 - Q = ~J drdx [Kj(r) (47r|(r - x)!)" 1 £*(x) + £»(x) (47r|(r - x)))" 1 Kj(r)] . (3.8) 

Whether our observations about the behavior of quark packets coupled to each other by 
Hq_q also apply to quark packets coupled to gluonic matter directly by H g _g, depends on 
the distribution of glue — whether, when the integration in Eq. (|3.8|) ranges over regions in 
which x^r, the expansion in Eq. ( |3.7| ) is valid. We have to defer more detailed discussion of 
this question until the dynamical equations that determine the distribution of glue in this 
gauge-invariant formulation of the theory have been solved. 
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Another point of interest about Eq. (|3.6| ) is the fact that the spatial dependence of 
the long-range nonlocal interaction described by Hq_q is largely determined by jF 6a (r,x), 
which represents the expectation value jF ba (r,x) = {g\T (r, x)|gr). And \g) implements 
the Gauss's law given in Eq. (3.5), in which the Gauss's law operator includes only the 
gluon color charge density. We observe, therefore, that the gluonic medium has a far more 
important role in determining the long-range behavior of the nonlocal interaction Hq_q 
than the quarks, which are acted on by Hq^q (and by H 9 ^q), but which have no role in 
transmitting these interactions. This result can therefore connect our work to other models 
for understanding color confinement. If it turns out that quarks, and not only gluons, are 
essential in transmitting confining forces from one group of quarks to another, then Hq_q 
and Hg_Q cannot be the only mechanisms for quark confinement. However, if other methods 
of analysis corroborate that the gluonic medium has the primary role in effecting quark 
confinement, than Hq_q and H g _q become more interesting candidates as descriptions of 
confining forces. 



IV. NON-PERTURBATIVE APPROACH TO THE NONLOCAL INTERACTION 

IN AN SU(2) MODEL OF QCD 

In order to pursue the analysis undertaken in Section [TIT] , and to demonstrate a procedure 
for evaluating Hq_q nonperturbatively, we will make use of Yang-Mills theory — the SU(2) 
version of this model — for which the structure constants f 5ba are e 5ba . We will also ignore 
correlations among gauge- invariant gauge fields, and replace all A 5 Gli (r) in Eq. (|3~3| ) with 
the corresponding (A Gii (r)) = (g\ A Gii (r) \g). Furthermore, we will model (A Gii (r)), which 
are transverse fields in the adjoint representation of SU(2), as the manifestly transverse 
"hedgehog" configuration 

<^i<(r)> = e«V(r). (4.1) 

Although there is no reason to believe that the ansatz given in Eq. follows from the 

dynamical equations that determine A s Gii (r), it is a convenient choice for examining to what 
extent the structure of ^ fta (r,x) — as distinct from the precise form of A Gii (r) — enables 
us to nonperturbatively evaluate the infinite series given in Eq. ( |3.3j ). This simplified SU(2) 
model can also help us to identify the features of the gauge-invariant gauge field that might 
be significant for the confinement of quarks or color-bearing quark packets. 

When we substitute e"^ 7 for /S) , and use Eq. (|4.1| ) in Eq. ( |3.1|) , we find that the integral 
in Eq. (|3.1| ) that is linear in the gauge-invariant gauge field, can be expressed as 

(-9) e %)6a (4ii( r )}^ / 47r| ^ Xi| Jo(xi) = ige^ir)^ J (4.2) 

where L s represents the (^-component of the orbital angular momentum. In the case of the 
term quadratic in the gauge-invariant gauge field, we observe that 
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9 V«'-<.>e-<*>"<4,'*(r)>a, /sj^i <4,1(x 2 ))S, / i - j ^ f jS(x 1 ) = 

- 9 V(>>»". e «%)>(r)L iM / , / — , jS(x t ). (4.3) 

J 47r|r — x 2 | J 47r|x 2 — xi| 

L S( - 2 ~> can be shifted to the left of the x 2 integration, by noting that for any reasonably 
well-behaved ip(x 2 ), 



47r|r-x 2 | d(x 2 )j J \ d{x 2 )j 47r|r - x 2 |/ 

= -e^r—Jd^^x,) ( —-L— \ ^(x 2 ) = -il/>Jdx*<t>{x 2 ) f — j— ^(x 2 ), (4. 

and that the identical procedure can be carried out on every term in the series given in 
Eq. ( p.ip , to yield an expression in which all orbital angular momentum operators are on 
the extreme left-hand-side of the expression, and are functions of re and d/dry. We thus 
obtain 

oo oo 

KoM = E ^t) L U^U v ) = -io(r) + E i B «$4o*J B) (r) (4.5) 

n=0 ra=l 

where L^L = 112=1 an d is given by the SU(2) version of Eq. fl2.7|) (with the convention 



that eg 1 = e 5fea , eg } a = -5 ba , and Lf 0) = 1), and 



*(„,(r) = #(r) / - , dy " ) rf fo) / - . ^ . »(»>)■ ■ ■ / ^pp^ . . x 

' 4yr|r — y (1) | J 47r|y ( i) - y (2) | J 4vr|y^_i) - y w | 

dyi "- 1> r^^-D) / , ,., dy(, ' ) ., tf(y w ) , (4.6) 



47r|y(n-2) - y(n-i)| J 4tt |y (n _i) - y (n) 

with the convention that $" ^(r) = Jo(r). 

Because of the simplicity of the SU(2) structure constants — the Kronecker delta is 
the only symmetric structure constant in SU(2) — it is possible to significantly simplify 
i n e s ^L s ^ n y It was previously pointed out that the SU(2) chain of structure constants e^ 7 
can be represented as 



e (*n) 7 = ( _1 ) 2 1( W]«[2] ^«[3]a[4] " ' ' 5 a [n-3]a[n-2] ^ ^ e^ 7 (4.7) 

and 

e 1n) = (~ l)~^a[l]a[2] ^a[3]a[4] ' " " #a[n-2]a[n-l] e"^ 7 (4.8) 

for even and odd n respectively. M It is a trivial consequence of Eqs. (|4.7| ) and ( |4.8| ) that 

<5ba r <5 
(n)^(n) 



jnjbaj^s can k e represented as 



i n ef^Lf n) = Cfc = A n (L a L b + L h L a ) + B„5 afe + zC n e a&i L^' , (4.9) 

where A n , B„, and C n only depend on numerical constants and on the Casimir operator 
L 2 . From simple recursion relations, i.e. D n + D( n _i) = D( n _ 2 )L 2 where D n = A n + C n = 
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B( n+1 )/L 2 , and F„ + F( n+1 ) + D n + D( n _x)L 2 = where F n = A n — C„ , we can obtain the 
following explicit expressions for A n , B n , and C n . 



\n-2l 



An = Lil! V (n ' * )! L 2s for n>2, (4.10) 

2 ^ (n-2a-2)!(a + l)! 1 ; 



— 

1 ~ ; 1 n-s-2)! 



K = {-lt~ l) E 7 V 9 9 / 2(S+1) for n>2, (4.11) 
(n-2s-2)!(s)! 



and 



rn-ll 



2 ^ (n-2s- l)!(s + l)! 

where [(n — 2)/2] = (n — 2)/2 when n is even, and [(n — 2)/2] = (n — 3)/2 when n is odd. 
for n < 2, A = Ai = Bi = Co = 0, and B = Q = —1. Because the algebra of the elements 
(L a L b + L b L a ), S a b and ie a bjLP is closed under multiplication, it is easily shown — using the 
notation introduced in Eq. ( f4.9|) — that 

£(n+fc) = — £(n)£(fc) (4-13) 

where the repeated index s is summed, and where 

A n +fc = (2L 2 — A n Afc + (A n B,t + B n Afe) — - (A n C k + C n A k ) — -C n C k , 
Bn+fc = A n A fc L + (A„C fc + C n A k ) L + B n B fc + C n C k L , and 

C„ +fc = - (2L 2 - ~) A n A fe + (B„C fc + C n B k ) + - {A n C k + C n A k ) - ^C n C k . (4.14) 

Since there are only three independent coefficients in any Cfy - viz. A n , B n and C n — the 
following equations must always have a unique solution for the constants for i = 1— >3: 

- A = X1A1 + x 2 A 2 + X3A3 , 

— B = X1B1 + x 2 B 2 + X3B3 , and 

-C = X1C1 +x 2 C 2 + X3C3, (4.15) 

so that the linear combination 

C% + x 1 £^ ) + * 2 C\% + xs£fa = . (4.16) 

Moreover, Eq. (|4.13|) has the effect of generalizing the validity of Eqs. ( [4.15 ) and ( 4.16| ) so 
that they apply to coefficients A n+i , B n+i and C n+ j for any fixed n and i = 0—>-3. The effect 
of these identities and of the explicit values of the three xi is that 

£(n+3) + 2£(„ +2 ) + (l — £ 2 ) — L 2 ^^ = (4-17) 
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for any n. The Casimir operator, L 2 , is treated like any constant in this analysis. 



Eq. (|4.17f) enables us to obtain a differential equation for /Cg(r) whose features — in- 
cluding its order — reflect the SU(2) algebra as expressed in Eq. (fO|). We observe that, for 
n>l, 



! 1 frb 

0(r) (n) 



which, with the use of Eq. (|4.13|) , leads to 

_o 1 



<j)(r) 



K(r)+j b (r) = £&/Cg(r), 



(4.18) 



(4.19) 



(r) 



V 2 



^(r)+j b (r) 



>L (1) V 



0(r) 



j*(r) = -£&/C3(r) and (4.20) 



v 2 4^(v 2 1 



(r) y 0(r) 



^(r)+j 6 (r) 



-4 S ) V2 ^o(r) f = ^3)^o(r) 



(4.21) 



We can use Eq. ( [4.17 ) to combine Eqs. ( 4.19|) - (fi.21 ) in such a manner that the right-hand 
sides vanish, to generate a sixth-order equation for K.q(t), given below: 



v 2 — (v 2 - 1 *" b 



*3(r) 



_2V 2 — IV 2 - 1 ^ 



^(r)]+L 2 ^(r) + 



+ [25 ba -C b ^)V 



V 2 



:Jo(r) 



V 2 



V 2 



V 2 



V 2 



yJo(r) + 
rJ'o(r) 



(4.22) 



Although Eqs. ( |4.5|) and (|4.6| ) represent /Cg(r) as an infinite series, which we might expect 
to have to evaluate perturbatively, Eq. ( |4.22p is a differential equation that JCq(t) obeys as 
a whole. The derivation of Eq. ( |4.22|) therefore has enabled us to bypass the perturbative 
representation of /Cg(r). 

We will make the further ansatz that <fi(r) is a complex constant, k. Then, we can write 
Eq. (|4.22|) in the form 



V " 2^ + (1 - L 2 )V 2 + kL 2 ) JC b (r) = s b (r) 



K 2 



(4.23) 



where s (r) is the source term 



s r 



1 - L 2 ) 6 ba + (2£ft + 4 a ))] V 2 + (25 ba - £*) — - — | j b (r). (4.24) 
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In Eq. ( |4.24| ), s b (r) has been expanded in inverse powers of k, and the leading term has 
been kept ^-independent. Eq. ( [4.23|) is a differential equation in which /Cg(r) is related to 
a source. The spatial dependence of /Cg(r) can be determined by finding a Green function, 
G(r, x), for which 



'V 6 V 4 \ 
-V - 2— + (1 - L 2 ) V 2 + kL 2 G(r, x) = -5(r - x). 



(4.25) 



Because the differential operator in Eq. ( f4.25p is explicitly dependent on the orbital angular 
momentum L 2 , it is desirable to expand G(r, x) in terms of partial waves. We express G(r, x) 
as 



Ye 



m \ "Xi fx J - 1 z,m\ w ri fr 



and find that ge(r,x) obeys 



(4.26) 



V r gt(r,x) = ~^5(r - x) 



where T> r is the differential operator 

3 c / 



V - 



+ 1 



- - V - 



£{e+iy 



K 



/V 



+ [!-£(£+!)] V r 2 - 



where V 2 = d 2 + {2/r)d r . We expand g^(r, x) as 



2 f°° 

ge{r,x) = - ge(a)ji(ar)j e (ax)a 2 da 



7T JO 



and adjust gi(a) so that 



2 r°° 1 

V r g e (r,x) = / ji(ar)j e (ax)a 2 da = — -6(r - x), 

ix Jo r z 



which leads to 



K 



a 



' ot A a 2 



- + 1 7^ + --^ + !) 



K 



(4.27) 



+ nl (£+1) 
(4.28) 

(4.29) 
(4.30) 
(4.31) 



Since k is a complex constant, we can parameterize it as k = |re|exp[z/3] and, for I > 0, 
represent gt(r, x) as 



K 2 r°° 
(r, x ) = — 

71 J-oo 



ji(ar)j((ax)a 2 da 



a 



(a - iyJ\K\e^/ 2 ^ (a + i^f\K\e^/ 2 ^ 
- ae^' 2 ) (a + ae^ 2 ) (a - ih^' 2 ) (a + ibe^' 2 



-i 



(4.32) 



where 
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\ 




+ 1) + 1 



and b 



\ 




(4.33) 



For £ = 0, the corresponding expression is 



9o{r,x) 



K 



nrx 



sin ( ar ) sin (ax) da 



00 a 2 



a 



K\e 



i/3/2 



) ( 



a + i\ \k e 



,1/3/2 



(4.34) 



The integrals in Eqs. ( 4.32 ) and ( 4.34j) are simple to evaluate, but only some of the features 
of the resulting expressions are relevant to this discussion. For I > 0, the dominant behavior 
of all the ge(r, x) as r and x increase in size is to decay exponentially in r and x. For £ = 
the poles at a = ii^/j^le^ 2 result in contributions that similarly decay exponentially in r 
and x unless e 4/3 / 2 = ±i; but the doubly degenerate pole at a = produces the contribution 



[9o(r,x)) 



(q=0) 



2rx 



((r + x) — \r — x\) . 



(4.35) 



The discussion presented in this section leads us to make the following observations: 
The series representation, in Eq. ([4.5[), of /Cq(t) — a quantity that encapsulates the nonlocal 
interactions between quark color-charge densities Jo(r) with each other and with gluonic 
color-charge — together with the SU(2) identities given in Eq. ( |4.17|) , lead to a differential 
equation for JCq(t) with source terms that are functionals of jg(r). The derivation of this 



differential equation — Eq. ( 4.22j) — enables us to eliminate the need for a term-by-term 



iterative expansion of ICl(r). We have arrived at Eq. ( [4.22| ) by replacing the SU(3) structure 
constants that apply to QCD with their corresponding SU(2) equivalents, and we have 
simplified Eq. ( ^4.22| ) to the form given in Eq. ( |4.23| ) by imposing an ad hoc ansatz that 
fixes the functional dependence of the gauge-invariant gauge field on spatial variables and 
SU(2) indices as shown in Eq. ([4.1[). In spite of the special assumptions that apply to this 
nonperturbative evaluation of /Co(r), it can serve as a useful model for a similar approach 
applicable to a more realistic treatment of QCD with SU(3) structure constants and with 
gauge-invariant gauge fields that reflect more of the the dynamics of this theory. 



The solutions we obtained for our simplified SU(2) model — Eqs. ( |4.29| )-( [4.35| ) — indicate 
that our simplifying assumptions lead to a form for the gauge-invariant gauge field that does 
not make Hq^q a confining nonlocal interaction. The fact that Eq. ( |4.23| ) is a sixth-order 



equation might have led us to anticipate that J 7 (r, x) would be a confining interaction. 
Even fourth-order equations can lead to Green functions with linear potentials that confine. 
2"7f But, in this present case, the SU(2) structure constants lead to the pole structure for 
the £-th partial wave shown in Eq. (|4.32|) , with ^Jn£(£ + 1) acting like a mass term in the 
differential equation that defines the Green function. Green functions that confine with 
linear potentials, or even with more rapidly rising ones, require higher order degenerate 
poles. Moreover, if the Green function's exponential decrease with distance is to be avoided, 
the degenerate poles must be on the real-a axis; and, if oscillatory behavior of the Green 
function is also to be avoided, the degenerate poles must be at the origin in the a plane. In 
our model, only g (r,x) has a degenerate pole at the origin, and the pole is not degenerate 
to a sufficiently high order to support confinement. We could expect that the richer algebra 
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of the SU(3) structure constants would lead to a higher order differential equation, and that 
the SU(3) Green functions therefore could lead to confinement — perhaps through a more 
degenerate pole structure. But in order to determine whether that is the case, it will be 
necessary to find an expression for the gauge-invariant gauge field that adheres more closely 
to the dynamics of the theory, and is not based on an ad hoc ansatz. 

To illustrate this discussion with a specific example, we consider a hypothetical case 
in which both L 2 and 1 — L 2 in Eq. ( |4.23[ ) vanish — an obvious impossibility, since the 
eigenvalues of L 2 are quantized to a set of possible eigenvalues that forbid this. Nevertheless, 
the consequences of this assumption serve a useful illustrative purpose. For this hypothetical 
case, we obtain the following expression for g (r,x), the £ = component of the Green 
function: 

k 2 f°° sin(o;r)sin(aa;)(ia 

= o*(o» + ao - (4 - 36) 

which includes a contribution from the quadruply degenerate pole at a = 0, 

[9o{r, x)] (a=0) = ((r + x) 3 - \r - x\ 3 ) . (4.37) 

The spatial dependence of [go(r, a:)]( a =o) is consistent with a confining potential for color- 
bearing quark packets. Although the dynamics of the model we are exploring does not 
lead to Eq. ( |4.37| ), the result may nevertheless serve a useful illustrative purpose. The 
model we are investigating in this section, which includes SU(2) structure constants and the 
"hedgehog" representation of the gauge-invariant gauge field, with the spatial function </>(r) 
represented as the constant 4>{r) = k, is itself a toy model used to represent the less tractable 
theory that has SU(3) structure constants and a gauge-invariant gauge field obtained from 
the dynamical equations presented in Refs. It is relevant to inquire what changes can 

occur in the behavior of /Co(r), and therefore also of the coupling term for quark color- 
charge densities, JF 6a (r,x), when changes are made in the equations that determine the 
Green function G(r, x), which might well be duplicated in the full SU(3) version of QCD. 

The model proposed in this section, consisting of SU(2) structure constants and the 
hedgehog ansatz for the gauge-invariant gauge field, does not provide us with much in- 
formation about Hg_Q — the direct coupling between quark and glue color. When the 
gauge-invariant gauge field appears as part of jF 6a (r,x), it is evaluated in expressions that 
do not contain the operator II£(r) conjugate to the gauge-dependent gauge field. Ignoring 
field correlations, and replacing Aq H (t) with (A S Q H (r)) in jF ba (r, x), therefore is a legitimate 
and useful approximation. The model is not, however, as applicable to the representation 
of Kj(r) = gf dae Tr [^^(rjf y c (r)f] A^^U^r), in which A 5 GH (r) and nj(r) appear to- 
gether. Our model does not include an expression for II^(r) that would be consistent with 
the hedgehog ansatz, or even respect the necessary commutation rules between the gauge 
field and its canonical adjoint. The effect of H 3 _q on color confinement therefore remains 
to be addressed until the dynamics of the gauge-invariant gauge field have been more fully 
explored. 
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V. DISCUSSION 



In this work, we have analyzed the nonlocal interaction that results when QCD is for- 
mulated in terms of gauge-invariant quark and gluon operator-valued fields. We have shown 
that this nonlocal interaction involves the quark color-charge density in a way that is roughly 
analogous to the role of the electric charge density in the Coulomb interaction in a gauge- 
invariant formulation of QED; but the functional form of this nonlocal interaction — the 
QCD analog of 1/ (47r[r — x|) in QED — is jF 6a (r, x), which is a nonlocal functional that 
depends not only on spatial variables, but also involves the gauge-invariant gauge field. 
jF ba (r, x) is an infinite series, in which the n-th order term contains the gauge-invariant 
gauge field to the n-th power. But the series has such a regular structure, that an explicit 
form for the n-th order term can easily be written, without requiring knowledge of the 
lower order terms in the series. A nonperturbative treatment of this nonlocal interaction 
term is therefore not nearly as inaccessible as would be the case for S-matrix elements in 
perturbative QCD. 

One feature of the nonlocal interaction between quark color-charge densities is that the 
monopole color charge is its leading term, and that higher order multipole interactions - 
Van der Waals type forces — succeed it in a series expansion. A color-charge monopole force 
as the leading term for quarks or color-bearing ensembles of quarks, and the consequence of 
that idea — that color singlets, the natural QCD analogs of electrically neutral ensembles 
of charges, would not feel that force — is not new. This idea — and the term "color 
neutrality" to designate it — have, for example, been suggested by Jain, Pire and Ralston. 
|26j These authors also have suggested that there is a connection between color transparency 



and the lesser importance of higher order Van der Waals-like multipoles for small-sized color 
singlet configurations. What is new in our work is that these features of a nonlocal QCD 
interaction are no longer a conjecture motivated by phenomenology only, but are inherent 
in the representation of the QCD Hamiltonian in terms of gauge-invariant quark and gluon 
fields. 



The gauge-invariant form of the QCD Hamiltonian contains the kinetic energy terms 
for the gauge-invariant quark and gluon fields, the nonlocal interaction discussed in this 
paper, and the additional interaction term Hj-A = —g J 4'K r ) a i^' l P( r ) y ^c\ j( r )^ r ? "which is 
the QCD analog of the QED interaction term — e / ip Jf (r)aii(j(r)AQ\ j (r)dr. In QED, this 
latter term describes the interaction of the electron current with the gauge-invariant gauge 
field (in QED, Aq\i (r) is just the transverse component of the gauge field). In QED, this 
interaction couples electrons to the two helicity modes of the photons — the propagating, 
observable quantized modes of the transverse electromagnetic vector potential. As is well 
known, the current appearing in this interaction has a v/c dependence, that makes it rela- 
tively unimportant in the low-energy regime, in which the Coulomb interaction is, by far, 
the most important interaction between charged particles. Since the corresponding term in 
the QCD Hamiltonian similarly involves a current density — the transverse color-current 
density in this case — it is reasonable to expect that its v/c dependence will also make it 
much less important than the nonlocal interaction between quark color-charge densities in 
the low-energy regime. The interaction term H g _q + Hq_q therefore is a very interesting 
candidate low-energy limit of of the QCD interaction — a nonlocal interaction term between 
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quark color-charge densities with each other and with gluonic color, that, in analogy with 
the Coulomb interaction in QED, would describe the most important features of QCD in 
the low-energy regime. 

It is important to make the following distinction between different parts of this work: 
On the one hand, there are relations like Eqs. ( |3.1| )-( ]3~3| ), which are exact. These relations 
are useful to the extent to which we have identified serviceable gauge-invariant fields that 
lead to a gauge-invariant Hamiltonian in which the nonlocal interactions that result from 
the implementation of gauge invariance are dominant in describing static quarks. These 
relations are not dependent on any approximations. They are direct and inevitable conse- 
quences of transforming the QCD Hamiltonian to a representation in which it is expressed 
in terms of the gauge-invariant fields we have constructed. On the other hand, there are the 
contents of Section [TV], which serve an important illustrative purpose, but are dependent on 
a number of simplifying assumptions. Our discussion, in Section |T|, of the implications of 
the nonlocal interaction for color confinement and for color transparency, are not dependent 
on the substitution of the SU(2) for the SU(3) algebra, or the ansatz that the gauge-invariant 
gauge fields are uncorrelated or that they have a hedgehog configuration, which we made 



in Section [TV]. This part of the discussion in Section |T| is dependent only on assumptions 
about the functional dependence of ^ fea (r,x) on the spatial variables r and x , and about 
the sizes of the quark packets coupled by this nonlocal interaction. We have developed 
this model far enough so that we can explore the consequences of our model, given these 
assumptions. We cannot, at this point, establish that these assumed quark configurations 
will necessarily arise in the low-energy regime. 



The differential equation, Eq. ( |4.22j ) and its special form Eq. ([4.23|) , and the Green 



function solution of Eq. (|4.23| ) given in Section [IV], depend on two simplifying assumptions 



- the SU(2) algebra and the hedgehog configuration of the gauge-invariant gauge field. 
These results are important for establishing a pattern for nonperturbative treatments of 
the nonlocal interaction in QCD and for demonstrating the role of the SU(N) algebra in 
generating higher order differential equations for the Green function that connects quarks 
to each other. But the specific solutions — Eqs. (}4.32|) - (|4.35|) — reported in Section III only 



apply to the "toy model" based on SU(2) structure constants and the hedgehog ansatz for 
the gauge-invariant gauge field. 



ACKNOWLEDGMENTS 



This research was supported by the Department of Energy under Grant No. DE-FG02- 
92ER40716.00. 



20 



REFERENCES 

[1] M. Belloni, L. Chen and K. Haller, Phys. Lett. B 373 (1996), 185. 

[2] L. Chen, M. Belloni and K. Haller, Phys. Rev. D 55 (1997), 2347. 

[3] M. Belloni, L. Chen and K. Haller, Phys. Lett. B 403 (1997), 316. 

[4] This material is reviewed in K. Haller and E. Lim-Lombridas, Found, of Phys. 24 (1994), 
217, and further references are cited therein. 

[5] K. Haller and E. Lim-Lombridas, Ann. Phys. (N.Y.) 246 (1996), 1; Erratum, ibid. 257 
(1997), 205; 

[6] P. E. Haagensen and K. Johnson, On the wavefunctional for two heavy color sources in 
Yang-Mills theory, MIT-CTP-2614, frep-th/ 97022031 

[7] M. Belloni, An Operator Approach to Gauss' Law in Yang-Mills Theory and Quantum 
Chromodynamics. University of Connecticut Ph.D. dissertation, 1997; (unpublished) 
Appendix F. 

P. A. M. Dirac, Can. J. Phys. 33 (1955), 650. 

[9] M. Lavelle and D. McMullan, Phys. Lett. B 329 (1994), 68. 

[10] E. Bagan, M. Lavelle, B. Fiol, N. Roy and D. McMullan, Constituent Quarks from 
QCD: Perturbation Theory and the Infrared, [hep-ph /960933C 



[11] J. Goldstone and R. Jackiw, Phys. Lett. B 74 (1978), 81. 
[12] R. Jackiw, Rev. Mod. Phys. 52 (1980), 661. 

[13] K. Haller, Phys. Rev. D 36 (1987), 1839; Phys. Lett. D B 251 (1990), 575. 

[14] G. C. Rossi and M. Testa, Nucl. Phys. B 163 (1980), 109; B176 (1980) 477; 
B237 (1984) 442. 

[15] K. A. James and P. V. Landshoff, Phys. Lett. B 251 (1990), 167. 

[16] D. Karabali and V. P. Nair, Nucl. Phys. B464 (1996), 135 

[17] K. Johnson, L. Lellouch, and J. Polonyi, Nucl. Phys. B 367 (1991), 675. 

[18] P. E. Haagensen and K. Johnson, Nucl. Phys. B 439 (1995), 597. 

[19] M. Lavelle and D. McMullan, Phys. Rept. 279 (1997), 1. 

[20] D. Stoll, Phys. Lett. B 336 (1994), 518; Phys. Lett. B 336 (1994), 524. 

[21] F. Lenz, H. Naus, and M. Thies, Ann. Phys. 317 (1994), 233. 

[22] F. Lenz, E. Moniz, and M. Thies, Ann. Phys. 242 (1995), 429. 

[23] T. D. Lee, "Particle Physics and Introduction to Field Theory", Chapter 18, Harwood 
(NY 1981). 



21 



[24] N. H. Christ and T. D. Lee, Phys. Rev. D22 (1980), 939. 

[25] G. A. Miller, Color Transparency - ■ Color Coherent Effects in Nuclear Physics. 
DOE/ER/41014-26-N97; |nucl-tli/9707040| . 

[26] P. Jain, B. Pire, and J. P. Ralston, Physics Reports 271 (1996), 67. 

[27] P. D. Mannheim and D. Kazanas, Gen. Rel. Grav. 26 (1994), 337. 



22 



